Reynolds number dependence of drag reduction by rod-like polymers 
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We present experimental and theoretical results addressing the Reynolds number (Re) dependence 
of drag reduction by sufficiently large concentrations of rod-like polymers in turbulent wall-bounded 
flows. It is shown that when Re is small the drag is enhanced. On the other hand when Re increases 
the drag is reduced and eventually the Maximal Drag Reduction (MDR) asymptote is attained. The 
theory is shown to be in excellent agreement with experiments, rationalizing and explaining all the 
universal and the non-universal aspects of drag reduction by rod-like polymers. 



I. INTRODUCTION 

The phenomenon of turbulent drag reduction describes 
the increase in the throughput in turbulent flows by 
adding tiny amounts of polymers. In this way for in- 
stance the amount of liquid that is transported in a pipe 
for a given pressure drop can be increased significantly. 
This implies a wide range of industrial applications, and 
consequently the effect, discovered in 1944 by Toms pj, 
has been studied intensively over the past decades (see 
e.g. Ref. 2]). The aim of this paper is to rationalize 
and understand the change from drag enhancement to 
drag reduction by rod-like polymers in turbulent wall- 
bounded flows as a function of the Reynolds number Re 
0, [J 01 ■ To this aim we present new experimental data 
and a theoretical analysis. 

A classical example of the phenomenon of interest is 
shown in Fig. ^ which pertains to pipe flow 4j. We 
denote the velocity field U(r,t) and its mean over time 
as V(y) where y is the distance from the wall. With the 
mean shear defined as S(y) = dV(y)/dy, the Fanning 
drag coefficient is defined as 

f^T w /(\pV*) , (1) 

where t w is the wall shear stress at the wall 




puS(y = 0) , 



(2) 



v, p and V are the kinematic viscosity, the fluid den- 
sity and the mean fluid throughput, respectively. Fig. Q] 
presents this quantity as a function of Re in the tra- 
ditional Prandtl-Karman coordinates l/y/J vs. Re^/7- 
The straight continuous line denoted by 'N' presents the 
Newtonian universal law. Data points below this line 
are indicative of a drag enhancement, i.e., an increase in 
the dissipation due to the addition of the polymer. Con- 
veresely, data points above the line correspond to drag 
reduction, which is always bound by the Maximal Drag 
Reduction (MDR) asymptote represented by the dashed 
line denoted by 'M'. This figure shows data for a rod-like 
polymer (a polyelectrolyte in aqueous solution at very 
low salt concentration) and shows how drag enhancement 
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FIG. 1: The drag in Prandtl-Karman coordinates for the wa- 
ter with the rod-likepolymer PAMH B1120 in 0.0001N NaCl 
in a pipe flow, see Q| for details. The symbols represent the 
concentrations in wppm (weight parts per million) as given in 
the figure. 



for low values of Re crosses over to drag reduction at large 
values of Re (3, 13 ■ One of the results of this paper is the 
theory presented below, that reproduces the phenomena 
shown in Fig. [l]in a satisfactory manner. 

It is interesting to note that with flexible polymers the 
situation is very different, and there is no drag enhance- 
ment at any value of Re. The reason for this distinction 
will be made clear below as well. 

In Sec. 2 we present new experimental results in which 
drag reduction is measured as a function of the Reynolds 
number, and for different concentrations of the rod-like 
polymer. We determine the rheological properties of 
these solutions simultaneously, so that the rheology of 
the solutions is known. In Sec. 3 we review the basic 
theoretical approach to drag reduction in general, and 
by rod-like polymers in particular. We discuss the effect 
of the polymers on the fluid flow in three cases: high Re 
flows, low Re flows and intermediate Re flows, respec- 
tively. In Sec. 4 we explain how to solve the equations 



and present the results and their comparison to the ex- 
periments. Discussions and the conclusions are presented 
in Sec. 5. 



II. EXPERIMENTAL 

A. Flow Geometry and Polymer 

In this section we describe the experimental results 
obtained using well-characterized rod-like polymer solu- 
tions. The turbulence cell has been previously described 
in detail by Cadot et al. 6]. A turbulent flow is gener- 
ated in a closed cylindrical cell (of volume three liters) 
between two counter-rotating disks spaced one disk di- 
ameter apart. The disks are driven by DC motors. The 
motors are configured to keep the disks rotating at a 
constant angular velocity f2, independently of the torque 
exerted by the turbulent fluid on the disks. Thus, all the 
experiments presented here are performed at constant an- 
gular velocity. All experiments were repeated with pure 
water and with the water seeded with the polymers. In 
both cases we define |(| the integral Reynolds number as 
Rc= pQR 2 /r) s , where p is the density and R is the ra- 
dius of the disks. To determine the energy dissipation in 
the turbulent fluid (from which we also infer the amount 
of drag reduction or drag enhancement), we take advan- 
tage of the fact that we are in a closed system. All the 
kinetic energy fed into the flow will eventually be dis- 
sipated viscously, leading to a temperature increase of 
the fluid: by measuring the temperature increase in time 
AT/ At we can also estimate the power dissipated by the 
turbulent flow. The temperature is measured with a Pt 
thermocouple probe with an accuracy of 0.01°C. From 
our temperature measurements, we can then deduce the 
percentage of drag reduction (DR) from: 



DR{%) ee 100 x 



(AT/At) poly 



(AT/At) water 



(3) 



For the polymer we use solutions of Xanthane, a stiff 
rodlike polymer with an average molecular weight of 
M w = 3.10 6 g.moP 1 . The rheology of the polymer solu- 
tions was studied on a Reologica Stress- Tech (cone-plate 
geometry) rheometer. The latter is equipped with a nor- 
mal force transducer, and has a large cone (55 mm) with 
a small angle (0.4°) in order to be able to detect small 
normal stress differences at high shear rates. Chang- 
ing the concentration allows us to modify the shear- 
dependent viscosity of the solutions. The apparent vis- 
cosity decreases with increasing shear rate, i.e., the fluid 
is shear thinning. A satisfactory description of most of 
the data for high shear rates can be obtained using the 
Ostwald-de Waele power law model pf for the viscosity 
ly. 7] = fciS^™ -1 ); the viscosity shows a power law depen- 
dence on the shear rate S. The validity of this model is 
limited to a certain range of shear rates, depending on 
the concentration C of the polymer solution; the range 




FIG. 2: Droplet detachment for a 300 wppm Xanthane solu- 
tion filmed with a rapid camera: there is 1ms between subse- 
quent images. 



and constants k± and n are reported in |8( ■ For the whole 
concentration range studied here (0 < C < 2000 wppm) , 
no measurable normal stresses were found in the range of 
shear rates studied (10 < S < 60008^); the uncertainty 
on the measurements is of the order of lOPa. 

The other important difference in the rheology of the 
rigid polymers, compared to solutions of flexible poly- 
mers, is the elongational viscosity tje, which quantifies 
the resistance to stretching of a fluid clement. If suffi- 
ciently important, t/e can be measured by looking at 
droplet fission with a rapid camera (Fig. The dy- 
namics of the thinning of the filament that connects the 
droplet to the orifice can be used to obtain the elonga- 
tional viscosity since the stretching of the filament cor- 
responds to a perfectly elongational flow. Performing 
this experiment, we find that, for the rigid polymer so- 
lution, rjE is so low that the dynamics of the filament 
is nearly indistinguishable from that of pure water for 
which rjE =3mPas [Tof . 



B. Results 

In Fig. 13 the percentage drag reduction DR(%) is 
plotted as a function of the Reynolds number for three 
different polymer concentrations. Also in our measure- 
ment geometry and with our rod-like polymer at low Re 
a drag enhancement is observed, which smoothly trans- 
forms into a drag reduction at high Re. In this paper we 
rationalize and explain the results presented in both Fig. 
^and Fig. [3| To this aim we turn now to a review of the 
theory of drag reduction. 
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A. High Re flows 
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FIG. 3: Percentage of drag reduction (or enhancement) vs. Re 
for three different polymer concentrations. Inset: percentage 
of drag reduction vs. Re when the actual viscosity (including 
the polymer) is used in the definition of Re) . 



III. REVIEW OF THE THEORY OF DRAG 
REDUCTION BY ROD-LIKE POLYMERS 

In the presence of a small concentration of rod-like 
polymers, the Navier-Stokes equations for the fluid veloc- 
ity U(r,t) gain an additional contribution to the stress 
tensor: 



dU 

— + U-VU = -Vp + v Q V 2 U + V • <t 
ot 



(4) 



V • U 







The extra stress tensor er is due to the interaction be- 
tween the polymers and the fluid. For rod-like polymers 
it takes the form [Tll | : 



(Tab = 6v p n a n b {ninjS i: j) , 



(5) 



where v p is the viscosity contributed by the polymer 
in the limit of zero shear, n = n(t, r) is a unit vec- 
tor describing the orientation of the polymer, and S%j = 
dUj/drj is velocity gradient. The equation of motion for 
the orientations of the polymers are approximated by the 
rigid-dumbbell model which reads 



DTZab 

Dt 



-6~f(1l ab - ^) , 



(6) 



where 1Z = nn and 7 is the Brownian rotational fre- 
quency, proportional to temperature. 



The phenomenon of drag reduction in highly turbu- 
lent channel flows can be understood on the basis of the 
balance equations of the mechanical momentum and tur- 
bulent energy For a channel flow we choose the co- 
ordinates x, y and z as the lengthwise, wall-normal and 
spanwise directions respectively. The length and width of 
the channel are usually taken much larger than the mid- 
channel height L, making the latter a natural re-scaling 
length for the introduction of dimensionless (similarity) 
variables, also known as "wall units" [jjj]. Denoting the 
pressure gradient p' = —dp/dx we define the friction 
Reynolds number Re T , the normalized distance from the 
wall y + and the normalized mean velocity V + (y + ) (which 
is in the x direction with a dependence on y only) by 



Re r 



p'L/v , y+ = yRc T /L , V+ = V/, 



(7) 



The balance equations are derived on the basis of the 
Reynolds decomposition 

Ui(r,t) = V(y)5 ix + Ui (r,t) , (8) 
Sij{r,t) = S(y)S ix 5 jy + s y (r,t) , (9) 

In addition to the mean shear S(y), we need to to in- 
troduce the mean turbulent kinetic energy K = (u 2 ) /2 
and the Reynolds stress W = —(u x u y ). The momentum 
balance equation is obtained by averaging Eq. Q and 
integrating with respect to y, ending up with the exact 
equation: 



u Q S + W+{a xy ) =p'{L-y) 



(10) 



The physical meaning of this equation is transparent: 
the momentum generated by the fixed pressure gradient 
p'(L — y) is transported to the wall by the momentum 
flux (^Reynolds stress) W, and dissipated by the New- 
tonian viscosity voS. The mean stress induced by the 
polymer {{<J xy )) is an additional viscous dissipation due 
to the polymers. 

In wall-units Eq. I|10|l can be written in the form: 

S+ + W+ + (a+ y ) = (1 - y+/Re r ) , (11) 

where S+ = v a S/{p'L), W+ = W/(p'L), and a± = 
aij/(p'L). When y + not too large, e.g., in the log-layer, 
we neglect the second term on the RHS, approximating 
the RHS as unity. 

The balance equation for the turbulent kinetic energy 
is calculated by taking the dot product of the fluctuation 
part of Eq. J3J with u: 

d 

WS = —{UyU 2 + UyP - diyUi) + VQ (Sjj' Sj,j ) + ((Tjj- Sy ) . 

Also this equation is exact. We simplify it by noting that 
the first term on the RHS involves the spatial flux of 
turbulent energy which is known to be negligible in the 
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log-layer. The second term on the RHS represents the 
Newtonian dissipation which is modelled (in wall-units) 
as IT3: 



IT 



3/2 



(13) 



where a and b are dimensionless coefficients of the or- 
der of unity and si = v Sij/(p' L) and if + = K/(p'L). 
The physical reasons of this model are follows: When 
y + is small, the dissipation of the turbulent energy is 
dominated by the viscosity (the firth term) . When y + is 
large, the an eddy with size y + and energy K + (y) loses 
its energy in the time-scale of y/y/K, which is the eddy 
turn-over time scale. We note that the summation sign 
in Ea. (|13|l is merely a pseudo sum. It only means that 
when y + is small (large), the first (second) term is more 
important. Writing in wall-unit, Eq l|12f) becomes: 



K- 



(K+) 



3/2 



(y 



ip 



(14) 



Finally, we quote the experimental fact that in the log- 
layer W + and K + are proportional to each other: 



K+c 2 = W + 



(15) 



Experimentally, it was found that c sa 0.5 in the Newto- 
nian case, and c w 0.25 in the MDR. 

To estimate of the polymer term on the LHS of Eq. 1|11|) 
we first observe that for large Re when the rod-like poly- 
mers are efficiently aligned by the flow, R xx as 1. There- 
fore from Eq. we conclude that 



(16) 



In wall units the momentum equation reads 

(1 + v+R vy )S+ + W+ = 1 - JL , (17) 

where v+ = v p /v Q . 

To estimate the polymer term in the RHS of Eq. Ijl2|l 
(the energy dissipation due to the polymers), we use the 
results of |lj| where it was shown that 



K 

{(TijSij) ~ VpRyy—^ ■ 



(18) 



Putting into Eq. (|12|) and writing in wall-unit, we have 

(19) 



K+ (K + ) 3 / 2 
W+S+ ~a 2 (l + v+Ryy)—+b [ - 



IP 



Equations l|17|) and (|19|) implies that the polymers change 
the properties of the flows by replacing the viscosity by 



v e s = 1 + v+R 



(20) 



This means that polymeric flows can be reasonably de- 
scribed by the changing the "effective viscosity" of the 



solution. In the fully developed turbulent flow, it was 
shown in [l4| that R yy depends on K + and S + : 



R K 

yv ~ (y+S+) 2 



(21) 



It was argued in ^4| that for large Re, K + grows linearly 
with y + and thus the viscosity profile is linear. Further- 
more, it was shown theoretically in ^5| that, if the effec- 
tive viscosity varies linearly with y + , i.e., 



^cff 



l + a(y + -y+) 



then c must satisfy the relation 



1 



(22) 



(23) 



Here <5 *=s 6 is the width of the viscous sub-layer in the 
Newtonian flows. 

For a given u eS , the equations l(T7|l. lfIS)l. (JTSJl and 
form a complete set of equations for variables S + , W + , 
K + and c. This model has been studied extensively for 
Newtonian flows ^| and for the flows with linear viscos- 
ity profile jl5| . The resulting velocity profiles agree rea- 
sonably with the experimental results. It was found that 
in the viscous sub-layer, the model does allows for a solu- 
tion that is turbulent, and therefore K + = W + = and 
S + = 1. This also agrees quite well with existing exper- 
imental data. Moreover, an edge solution for a = 1/12 
was observed [13 . If a > 1/12, the model no longer has a 
solution corresponding to a turbulent flow and therefore 
the flow must be laminar. The special case a — 1/12 was 
identified as the MDR of the polymeric drag reduction. 



B. Low Re flows 



According to Eq. I|20|) , the value of z/ g depends on v+ 
and Ryy The value of is determined by the polymer 
properties such as the number of monomers, their con- 
centration etc., and thus v£ should be considered as an 
external parameter in the equation. The value of R yy , 
on the other hand, depends on the properties of the flow. 
In the case of laminar flow with a constant shear rate, 
i.e., K + = W + = and S + =constant, it was shown 
theoretically in [llj that: 



Ryy — 



2 i/3 



De 



2/3 



(24) 



where the Deborah number De defined by De = S/j. 
Thus, the effective viscosity is reduced if S is increased, 
and therefore the rod-like polymers solution is a shear- 
thinning liquid. Naturally, The value of De changes with 
Re. To clarify this dependence we consider the momen- 
tum equations Ea. (|10|l at y = in the Newtonian case. 



UqS = p'L 



(25) 
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Usually in experiments system size and the working fluid 
remain the same. Therefore, Vq and L are constants and 
so Re T depends on p'L only. According to 10, Re T grows 
as y/p'L and therefore 



De 



^Re 2 
1 L 2 T 



by Eq.J23). Putting into Eci. (P4"|l . we have 



v cVi = 1 



A 



p Rcl/ 3 



where A = vq/^L 2 is a constant. 



(26) 



(27) 



where v p = f p /g, 5 + = S + g and if + = K + . We can find 
the a-dependence by comparing Eq. I|30|) to 122|) and 
identifying the width of the viscous sub-layer y^ with 
a/c{pt). This stems from the continuity of S + at the 
boundary of the viscous sub-layer. This means that 



a/c(a) = yt 



(31) 



Combining, Eqs (|2*2*|) . (|2*3^l and l(3*T]) . the relationship be- 
tween v e f{ and c is: 



a — cS ( , a 

^cff = H 7— ( V 

ad V c 



(32) 



If we note also that Eqs. (|2*2*|l and I3T)|| can be written as 



C. Intermediate Re flows 

In the case of intermediate Re, we need an interpola- 
tion between Eqs. I|21|l and l|24l) . To do this we note that 
when y + is small, the solution of Eqs. I|15l) . i|17[l and 
(|19fl result in W + = K + = in the viscous sub-layer. 
This implies that the flow cannot be highly turbulent in 
the viscous sub-layer. Thus, it is reasonable to employ 
Ea. (|24|l as long as y + is small. On the other hand, as the 
upper bound of y + is Re T , when y + is large, it automati- 
cally implies that Re r is large. The laminar contribution 
is therefore negligible as it varies inversely with Re T . The 
effective viscosity due to the polymer is dominated by the 
turbulent estimate, Eq. (|21|l . To connect these two re- 
gions we simply use the pseudo-sum: 



^eff = 1 




(y + s+) 



(28) 



where g = 1 + f+A/Rc^ 3 . One can see that the limits 
for both high and low Re r are satisfied. 



IV. SOLUTION AND RESULTS 



with 



K+ = A\S + y + f 



A 2 = jjjgg - 1 



(33) 



(34) 



Using then Eqs. (|15|1 and (|33|l . we can rewrite Eqs. 117|l 
and ljl9|l as two equations for the two variables v e g and 
S+: 



and 



hsS + + c 2 A 2 (S + y+) 2 = l , 



( ?S+=i> eS (^-f+bA^ 

y + 



Equation (|36|l implies 



S+ = 



v e ff a" 



(y+) 2 (c 2 -bA) 



(35) 
(36) 

(37) 



Substituting Eq. ^ into Eq. (33) gives an equation for 



~2 ( " \1(„1 



{ ^f^. bA) + c ^ 2 ^ =(c 2 -bA) 2 . (38) 



A. Setting up the equations 

Having the expression for v c g, we have to specify the 
value of c in Eq. (fTK|l to complete the equations. This 
variable naturally depends on a and g. The latter de- 
pendence, however, can be eliminated by rescaling the 
dimensional variables. Define 



Finally, we can solve Eqs. (|38|l and l|32|l to get iy e s{y + ) 
for different values of v p . Then we can obtain S + and 
K + using Eqs. I|37|) and l|33l) respectively. Finally, we 
reexpress the variables in wall unit by using Eq. (|29|l . 



B. Comparison of analysis with experiments 



In these units Eq. I|28|l is written as 

v ctt — 1 + v v = 

P (y+S+) 2 



V 



(29) 



(30) 



To compare our analytical results with the experi- 
ments, we first solve the model with Ea. (|28|l for Re T =590, 
and for A = 0, i.e. the high Re limit. This demon- 
strates that the MDR is reproduced within the model, as 
is shown in Fig. 01 

We observe that V + is increasing with a constant rate 
for small y + and then follows the log-law of the MDR. 
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FIG. 4: The velocity profile as predicted by the model with 
u+ = 10 4 and Re T = 590. 



FIG. 6: f~ 1/2 as the function of log 10 (Re/ 1/2 ) with A = 1 
with various values of v„. 
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FIG. 5: The velocity profile predicted by the model with 
Re T = 590 and A = 1 for various values of . 



The results for A = 1 are shown in Fig. |SJ We note 
that for f p = 100 the velocity profile first follows the 
viscous profile and then joins the MDR until it begins 
to cross back to the Newtonian plug. For f p =1000 the 
viscous layer is increased in size, but then the velocity 
profile becomes parallel to the MDR but with a smaller 
intercept, signaling less efneicient drag reduction. Once 
v p reaches the value of 10 000, the viscous layer becomes 
rather large, of the order of 20. The velocity profile sill 
succeeds to run parallel the MDR, but with a much re- 
duced intercept. 

To compare the prediction of the model to the experi- 
mental results shown in Fig. ^we need to relate Re T to 



Re as used in that figure. The relation is 

Re_L = Re T Vf , (39) 

where is the velocity at the center of the channel. 
Fig. [S]shows J" 1 / 2 as a function of log 10 (Rei/ 1 / 2 ) with 
various values of v p . One can see that for large v p , we 
obtain a drag enhancement when Re is small. This en- 
hancement decreases with increasing Re and eventually 
for large enough Re there is drag reduction. The maxi- 
mum amount of drag reduction, however, cannot exceed 
the MDR. For small i/ p , the drag enhancement is nearly 
unnoticeable. However, drag reduction in this case is 
smaller for large Re because the Newtonian plug occurred 
earlier. It should be noted that the curves in Fig. ex- 
hibit the "ladder" characteristics for different values of 
v p . This was identified as "Type B" drag reduction by 
Virk and is the fundamental characteristic of drag 
reduction by rod-like polymers. We see that our results 
agree reasonably well with the experiments presented in 
Prandtl-Karman coordinates. 

To compare the theory to the experimental results in 
DR(%) vs Re coordinates we use the experimentally mea- 
sured values of g and v p to predict the amount of drag 
reduction. As we mentioned before, the effective viscos- 
ity in laminar flow of XG solutions is well-approximated 
by Ostwald-dc Waclc model. This implies in the low-Re 
flow, g is given by (in MKS units): 



g = hS n ~ 



L 2 



(40) 



where the last equality is obtained by multiplying Eq. 
(|26[) by 7. The experimental value of v p for XG solution 
is given by (in MKS units): 

v p = 0.011147 C 1A22 v , (41) 
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FIG. 7: f~ 1/2 as the function of log 10 (Re/ 1/2 ) with A = 1 
with various values of v v . 

where C is the concentration of XG in wppm. The value 
of / can then be calculated by: (i) obtaining the experi- 
mental values of C, fci and n for Eqs. I|4l)|l and (|4ip. (ii) 
rescaling the variables S, K and W by Eg (|29fl and finally, 
(iii) solving Eqs. (|38|) and l|32|l to get / as a function of 
Re^. It should be noted that in the present experiments 
the definition of Reynolds number is in general different 
from our definition in 139|) . For the sake of comparison 
of the theory with the experiments we assume that two 
Reynolds numbers are proportional to each other: 

Re L = 12Re cxp . (42) 

Fig. |7| shows the comparison of the percentage of drag 
reduction (enhancement) between the theoretical predic- 
tions and the experimental results. The two data sets 
shown pertain to c = 250wppm (fej = 11.04, n — 0.727) 



and c = 500wppm (ki = 11.04, n = 0.727). The agree- 
ment between theory and experiment is very satisfactory. 



V. DISCUSSION AND CONCLUSIONS 

In summary, we studied the Re dependence of the drag 
reduction by rod-like polymers. We showed that when 
the value of Re is small, the drag is enhanced due to the 
homogenous increase of the effective viscosity. When Re 
is sufficiently large, due to the turbulent activity, the ef- 
fective viscosity varies as a function of the distance to the 
wall. As the result, the amount of drag is reduced. We 
use a simple interpolation between low Re and the high 
Re flows to account for this. The resulting theoretical 
results agree semi-quantitatively with the experiments 
both in Prandtl-Karman coordinates and in DR(%) vs 
Re coordinates. 

It should be noted that only rod-like polymers exhibit 
drag enhancement for low Re. For flexible polymers, the 
drag is the same as that of the Newtonian flow for low Re, 
and onl y af ter a critical value of Re, the drag reduction 
sets in [ljl Hal- This difference in behavior from rod- 
like polymers is because the flexible polymers are coiled 
when Re is small. They do not affect the flow unless 
Re increases enough to allow the shear to develop to af- 
fect the coil-stretch transition in the flexible polymers. 
In contrast, rod- like polymers are always extended and 
therefore they can affect the flow for all Re. 
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